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ABSTRACT 

We reconsider the issue of proving large scale spatial homogeneity of the universe, 
given isotropic observations about us and the possibility of source evolution both in 
numbers and luminosities. Two theorems make precise the freedom available in con¬ 
structing cosmological models that will fit the observations. They make quite clear that 
homogeneity cannot be proven without either a fully determinate theory of source evo¬ 
lution, or availability of distance measures that are independent of source evolution. 
We contrast this goal with the standard approach that assumes spatial homogeneity 
a priori, and determines source evolution functions on the basis of this assumption. 
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1 INTRODUCTION 

Ever since the earliest cosmological models, the Einstein and 
de Sitter models, we have been trying to fit observations 
to the Friedmann-Lemaitre-Robertson-Walker (FLRW) spa¬ 
tially homogeneous and isotropic family of models. The suc¬ 
cesses of reproducing a Hubble redshift-distance law, calcu¬ 
lating the correct cosmic helium & deuterium abundances, 
and the prediction of a cosmic microwave background radi¬ 
ation (CMBR), have convinced us of its validity as a bulk 
description of the universe. However, proving that the ge¬ 
ometry of the universe is FLRW on the largest scales is not 
easy. In fact, the history of observational cosmology shows 
that each time improved instruments permit deeper surveys, 
the new data soon reveals inhomogeneities on the new scale. 

The best evidence for homogeneity comes from limits on 
anisotropy of both galaxy counts and the CMBR, obtained 
in each case by comparison of observations in different di¬ 
rections. However this is strictly speaking only evidence for 
isotropy about the earth; homogeneity follows only if we in¬ 
troduce a Copernican principle, either for galaxies or for the 
Cosmic Background Radiation (^; |^. Without this assump¬ 
tion, the models indicated are isotropic about us, but allow 
a spatial variation of the geometry and matter content that 
is spherically symmetric about our position.^ The Coperni¬ 


can principle is not really in dispute on a sub-horizon scalejTj 
but could be incorrect on a super-horizon scale if we accept 
theories such as chaotic inflation ((^; see also (P)). There¬ 
fore one would like to actually prove homogeneity for the 
observable region of the universe, rather than assuming it 
on principle, which is essentially what happens in the usual 
reproach. Similar issues have been discussed by Goodman 

There are several problems with demonstrating homo¬ 
geneity from observed data. The deeper observations are not 
only fainter and redshifted, they are also affected by proper 
motion, reddening and absorption due to interstellar mat¬ 
ter. These contribute to selection effects which are tricky to 
compensate for. But the main problem at large distances is 
the evolution of sources, since deeper observations are re¬ 
ceived from earlier cosmic epochs. Evolution can take place 
in source colours, luminosities, and sizes; at high redshift 
it can affect the type of source as well as their numbers. 
However in this paper, for simplicity, we shall only consider 
bolometric observations of one type of source. 

How does the number and brightness of the observable 
sources relate to the local density at different times? Recent 
evidence for a sharp fall off of the space density of quasars 
above a redshift of z = 5 ^ could be taken as ev- 


* In a universe that is isotropic but inhomogeneous, there are 
anthropic reasons why we might be near the centre, as argued in 


t Although one can construe the current uncertainties in the val¬ 
ues of cosmological parameters such as Hq and O as being evi¬ 
dence for different values on different scales — i.e. inhomogeneity, 
we are not claiming this here. 
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idence of inhomogeneity, though most attribute it to source 
evolution. The large population of faint blue objects found 
by sensitive optical surveys is thought to be young starform¬ 
ing galaxies at high redshifts, and therefore constitute evi¬ 
dence for evolution (|^; @) — if we assume the universe 

has a FLRW geometry. Without evolution, these observa¬ 
tions are inconsistent with that geometry. One difficulty is 
that redshifts of faint objects are scarce and difficult to ob¬ 
tain. A redshift of z ~ 2 was deduced 0 by combining 
number counts versus magnitudes in 3 colours with galaxy 
evolution models and cosmological models, and comparing 
with the few measured redshifts available. Again this re¬ 
quired the assumption of homogeneity.^ Indeed, deducing 
the effects of source evolution by comparing observations 
with predictions in a FLRW model is a standard technique. 
Similarly studies of a luminosity-size relation also assume 
a FLRW model — recent examples are ( jl^ ; However 
this cannot lead to certainty (|^). The new discovery of a ra¬ 
dio galaxy which is apparently an immature giant elliptical 
galaxy at z — 4.41 (p^), provides a more striking exam¬ 
ple of probable source evolution. The problem is to show 
that this is not rather evidence of spatial inhomogeneity, 
manifested in a change of the evolutionary history or the 
nature of objects observed at larger spatial distances from 
us. Also, claims of a periodicity on top of the Hubble law in 
the redshift-distance relation ( p^ ) indicate significant devi¬ 
ations from standard FLRW observational relations, which 
could be due to spatial inhomogeneities (e.g. (p^) or to a 
temporal variation in the cosmic expansion rate. 

If suitably smoothed observations are isotropic, the 
principal observations of discrete sources one can hope to 
make are the number counts n{z) of sources as a function 
of “distance”, conveniently taken as given by cosmological 
redshift z, and the magnitudes and angular diameters of 
sources, also as a function of redshift. If the assumed linear 
size and absolute luminosities of the sources are correct, the 
latter two give the luminosity and area distances R{z), which 
should be equal. This is rather fortunate since in practice it 
is often difficult to separate the two measures — one has to 
define an edge of a galaxy image in order to measure its lumi¬ 
nosity, and conversely, one often dehnes the edge relative to 
the central brightness; and both measures are signihcant in 
determining selection effects (piP. In any case, at the largest 
distances angular diameters cannot be measured, and it is 
the luminosity distance that is used. 

We consider two types of source evolution; absolute lu¬ 
minosity L{z), and mass per source m{z), i.e. total density 
over source number density, which represents evolution in 
source number via sources turning on, galaxy mergers etc. 
Since source evolution is likely to be determined as a func¬ 
tion of age r, these functions could usefully be expressed 
as L{t{z)) and m(T(z)); however it is analytically easier to 
solve the observational equations if they are considered as 
functions of the observable 2 . This also helps to emphasise 
that if large-scale inhomogeneity were in fact to occur, the 


i It turns out many of the sample are not at high redshift, but 
there is still a “high redshift tail” to the sample (|l5|).) 


age of the universe would vary with spatial position and so 
becomes difficult to handle.|^ 

Earlier work (section 15.3 of (p^), section 7 of (p^) 
showed that if the observational relations are isotropic and 
of the FLRW form, then the universe is indeed homogeneous, 
provided we can assume the matter stress tensor is that of 
pressure-free matter .p] However that analysis did not fully 
consider the effects of source evolution. 

In this paper we show that any given isotropic set of ob¬ 
servations n{z) & R{z), together with any given evolution 
functions L{z) and m{z), can be fitted by a spherically sym¬ 
metric dust cosmology — a Lemaitre-Tolman-Bondi (LTB) 
model — in which observations are spherically symmetric 
about us because we are located near the central world¬ 
ling Thus we show that any spherically symmetric obser¬ 
vations we may eventually make can be accommodated by 
appropriate inhomogeneities in a LTB model — irrespective 
of what source evolution may occurp^ Conversely we show 
that, given any spherically symmetric geometry and any set 
of observations, we can find evolution functions that will 
make the model compatible with the observations. 

The purpose is to demonstrate explicity — develop¬ 
ing the ideas in (^) — that the relationship between the 
large scale isotropy of observations and large scale cosmic 
homogeneity is weaker than is commonly assumed. Indeed, 
apart from any other problems, we can’t have a good demon¬ 
stration of homogeneity without observational tests of our 
source evolution theories that are independent of cosmo¬ 
logical model, or distance measures that are not influenced 
by source evolution. This emphasises the conclusion that if 
the demonstration of homogeneity depends on knowing the 
source evolution, and validation of source evolution theo¬ 
ries depends on knowing the cosmological model is homo¬ 
geneous, then neither is proved. Indeed if we do not make 
the FLRW assumption, our results can be used to determine 
the degree of inhomogeneity from the observations and any 
given source evolution functions. If we do make the FLRW 
assumption, they can be used to determine the source evo¬ 
lution functions required to make the observations compat¬ 
ible with that model. The latter is the way theory is usually 
run. The point of our paper is to emphasise that there are 
other options, and so such source evolution results should 
be viewed with caution. 


2 THE LTB MODEL AND ITS NULL CONE 

We here outline the metric and our notation and null cone 
solution; for more details in this notation see (^sl) . 

The general spherically symmetric metric for an irro- 
tational dust matter source in synchronous comoving coor- 

§ Since inhomogeneity has been introduced, it is even possible 
the evolution functions are also position dependent. 

^ The result should also obtain for barotropic perfect fluids, but 
not necessarily for imperfect fluids. 

II Despite claims to the contrary in the literature, this is a per¬ 
fectly possible situation (^. 

** This is within the spirit of the programme of observational 
cosmology set out by Kristian and Sachs in their fund ame ntal 
paper (|24[) and developed further by Ellis, Stoeger, et al (|22|). 
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dinates is the Lemaitre-Tolman-Bondi (LTB) ( [^ p 8 | ) 

metric 

ds^ =-dV + + R^{t,r)dQ,'^, ( 1 ) 

where R'{t,r) = dR{t,r)/dr, and dQ^ = d6^ + sin^ Odtj)^. 
The function R = R{t, r) is the areal radius, since the proper 
area of a sphere of coordinate radius r on a time slice of 
constant t is AttR^. Solving the Einstein field equations gives 
a generalised ‘Friedmann’ equation for R{t,r), 


Wire 

and an expression for the density 

. / X M'(r) 

4np{t,r) = ^ 




(3) 


R^{t, r)R'{t, r) 

Eq (^) can be solved in terms of a parameter 77 = ri{t, r): 

(4, 


wherep] 
Sir) = J 


S{r) 


2E{r), 

1 , 

-2E{r), 


M{r) 


{ cosh v — 1, 

(1/2)W, 

1 — cos T], 


{ sinh T] — rj, 

(1/6)77®, when 
?7 - sin 77 , 


E > 0 
E = 0 
E <0 


(5) 


for hyperbolic, parabolic and elliptic solutions respectively. 

The LTB model is characterised by 3 arbitrary func¬ 
tions of coordinate radius r. E = E(r) > —1/2 has a geo¬ 
metric role, determining the local ‘embedding angle’ of spa¬ 
tial slices, and also a dynamic role, determining the local 
energy per unit mass of the dust particles, and hence the 
type of evolution of i?. M = M{r) is the effective gravita¬ 
tional mass with comoving radius r. ts = tsV) is the local 
time at which R = Q, i.e. the local time of the big bang — 
we have a non-simultaneous bang surface. Specification of 
these three arbitrary functions — M{r), E{r) and ts( 7 ') — 
fully determines the model, and whilst all have some type 
of physical or geometric interpretation, they admit a free¬ 
dom to choose the radial coordinate, leaving two physically 
meaningful choices, e.g. r = r{M), E = E{M), ts = ts(M). 


2.1 The Observer’s Null Cone 

We now take R and p as given on the observer’s past null 
cone, and we wish to express the 3 arbitrary LTB functions 
in terms of them, so characterising the LTB model that fits 
the observations. 

We generalise the gauge choice used in ( p^ ) to the case 
where the spatial sections are in general non-flat, i.e. all 
values of E. Human observations of the sky are essentially 
a single event on cosmological scales, so we only need to be 
able to locate a single null cone; we don’t need a general 

tt Strictly speaking, the hyperbolic, parabolic and elliptic so¬ 
lutions obtain when RE/M > 0, = 0, & < 0 respectively, since 
i? = 0 at a spherical origin in both hyperbolic and elliptic models. 


solution. On radial null geodesics, ds^ = 0 = dO^ = dtj)^-, 
so from (|^) if the past null cone of the observation event 
{t = to, r = 0 ) is given by t = t{r), then t satisfies 


R\t{r),r)^_ R' 


( 6 ) 


VTT2E VI + 2E 

We will denote a quantity evaluated on the observer’s null 
cone, t = t{r), by a"; for example R{t(r),r) = R. Now if we 
choose r so that, on the past light cone of {to,r), 


R'{t{r),r) 


R' 


= 1 , 


(7) 


VVT2E VVT2E 

then the incoming radial null geodesics are given by 

t(f) = to-r. (8) 

and the 


With our coordinate choice ([^, the density 
Friedmann equation (^) become 

M' 


A-KpR = 


VT+2E 


R = ± 


2M 

R 


+ 2E. 


(9) 


( 10 ) 


The gauge equation is found from the total derivative of R 
on the null cone 


dr dr 

and with (H) and (|l 0 [) substituted, it follows that 


dR 


- - VTT2E = -R= - ± J — +2E . 


dr 


2M 


R 


( 11 ) 


( 12 ) 


When we solve this for 2E(r) by squaring both sides and 
rearranging, we get 

\ 2 


1 + 275= i 


f \ +1 


/l^l . (13) 


This expression will tell us under what circumstances (or for 
which regions) the spatial sections are hyperbolic l-|- 2 i? > 1 , 
parabolic 1 -I- 2E = 1 or elliptic 1 -|- 2E < 1, based on data 
obtained from the null cone. We now use the expression 
for the density on the null cone to find a linear, first order 
differential equation for M(r). Eliminating 1 -|- 2E between 
(@) and (^, we get 


dM 

dr 


( 




4TTpR 

dr 


( 


M = 


V 


2ttpR'^ 

dr 


+1 


By evaluating (|^ and (^ on the null cone we find 


cosh?) — 1 , 
(l/2)77^ 

1 — cos 77 , 



(14) 


(15) 


{ sinh fj — fj, 
(1/6)77®. 

77 -sin 77 , 


when 



(16) 
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and 

r(r) = i{r) — tB{r) = to — r — is(r) 


(17) 


can be interpreted as proper time from the bang surface to 
the past null cone along the particle world lines. Thus, with 
M given by and then E by ([l3|), we can solve for fj from 


y SR 


and (^, T(r) from 


M ^ 

£ 3/2 ? 


(18) 


(19) 


with (^^ again, and hence tsir) from ( 0 )- 


2.2 Origin Conditions 

At the origin of spherical coordinates, r = 0, where R{t, 0) = 
0 and R{t, 0) = 0 for all t, we assume that the density is non¬ 
zero, that the type of time evolution (hyperbolic, parabolic 
or elliptic) is not different from its immediate neighbour¬ 
hood, and that all functions are smooth — i.e. functions of 
r have zero first derivative there. Thus eq (^ tells us that 
RE/M and E^^'^ jM must be finite at r = 0, (^) shows us 
that i? —> 0 and hence M —> 0 and E ~ at r = 0. Eqs 
( pd] ) and (0) become 


dr 


= i?'|r=0 = 


dr 


= VI + 2E = 1, 


( 20 ) 


and thus i? ~ r to lowest order near r = 0. From (g) we find 

( 21 ) 


M' ~ Anpor^ , M ~ -Trpbr® 


and so 
E 


(^Vo) 


2/3 


( 22 ) 


We verify these origin conditions satisfy (|l4|) to order and 
( 0 ) trivially to order r°. 


2.3 Redshift-distance formula 

We use the fact that in the geometric optics limit, for two 
light rays emitted on the worldline at rem with time inter¬ 
val Stern = t'^(rem) — t~(rem) and observed on the central 
worldline with time interval Stab = t'^ (0) — t~ (0) 


l + z = 


Stob 

Stem 


(23) 


which means that, integrating along the light ray and ap¬ 
plying this to the log of (p^), the redshift is given by 


(24) 


ln(l-|-a) = / R'{t, r) / \/l + 2E dr 

Jo 

for the central observer at r = 0, receiving signals from an 
emitter at r = rem- 

We need to find the redshift 2 explicitly in terms of 
observables. We differentiate (0) with respect to r: 


R'R 


M' 


MR' 


E' 


(25) 


V1 + 2F R-s/1 -f 2E RVl + 2E ^/1 + 2E 
so when evaluated on the observer’s past null cone, we get 

M' 


R' 

V1 + 2F 


R. 


- 4^ + (Vi + 2£;)' 


_R^l + 2E R2 
Now, from (0, the derivative of \/l + 2E is given by 

d^R 


(26) 


(V1 + 2F)' = 


d^R 

dr'^ 


M' 


M 


dr^ 


R 


+ ^- V1 + 2E (27) 
dR \ R2 


dr 


dr 


so, after eliminating M' by substituting from equation 
it follows that 


R' 


V1 + 2E 


i ( 4-RpR - MpR^l -h 2E / — 

r\ 

d?R d^R -— , dk 


, . :rk\ , (dR 


(28) 


where we have used equation ( |l^ ) to provide the second 
equality. From (M) it now follows that 


dr 


;in(l + 2)1 = - 


d^R 

dr^ 


+ Anpk 


/ 


dR 

dr 


(29) 


which theoretically gives the redshift in terms of coordinate 
radius r, directly from k(r) and p(r), viz 


ln(l + z) = - 


d?R 

dr^ 


+ A-Kpk 


/ 


<m 

dr 


dr . 


(30) 


However, we will be given observations in terms of 2 , 
rather than the unobservable coordinate r. This will be ad¬ 
dressed in the next section. 


The incoming radial null geodesics are given by 


3 OBSERVABLES AND SOURCE EVOLUTION 


dt = —R'(t, r) / \/l + 2E dr , 


so for two successive light rays, — & -f, passing through two 
nearby comoving worldlines ta & rs = -f dr at times t//, 

^ ^B 

—R'{t'^,r) + R'{t~,r)^ 
y/l + 2E ^ 


d(St) = StB — St A = df^ — dt = 


Consequently 

dlnSt = --^[R'(t,r)] /V1 + 2E dr 


For simplicity we shall confine ourselves to one type of cos¬ 
mic source and only consider bolometric luminosities. We 
shall assume that the luminosity of each source can evolve 
with time, and that the number density of sources can also 
evolve. The former we write as an absolute bolometric lu¬ 
minosity L, and the latter we shall represent as an evolving 
mass per source, m, which gives the total local density when 
multiplied by the source number density. As mentioned, we 
assume isotropy about the earth (once our proper motion 
has been accounted for), and also that the post decoupling 
universe is well described by zero pressure matter — “dust”. 
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The particles of this dust are galaxies (or perhaps clusters 
of galaxies). This means we can use the simplest inhomo¬ 
geneous cosmology — the LTB metric, which is spherically 
symmetric and inhomogeneous in the radial direction only, 
and is written in comoving coordinates. 

The two source evolution functions are most naturally 
expressed as functions of local proper time since the big 
bang, L{t) and mij). However, in a LTB model the time 
of the bang may vary from point to point, so that the age 
of objects at redshift z is uncertain both because the bang 
time is uncertain and because the location of the null cone 
is uncertain. The proper time from bang to null cone will 
be a function of redshift, t[z), and the projections of the 
evolution functions on the null cone we will write as L and 
rh. Of course, t[z) is unknown until we have solved for the 
LTB model that fits the data. However, for the sake of sim¬ 
plicity, we will take L and m to be given as functions of z, 
to illustrate how the 3 quantities, cosmic evolution, cosmic 
spatial variation, and source evolution are mixed together in 
the luminosity and number count observations, i and n. A 
treatment dealing with evolution functions based on r would 
involve solving a much more complicated set of differential 
equations in parallel. 


3.1 Relating Observables to the LTB Model 

The area distance or equivalently the diameter distance is 
the true linear extent of the source over the measured an¬ 
gular size. This is by definition the same as the areal radius 
in the LTB model R, which multiplies the angular displace¬ 
ments to give proper distances tangentially. The projection 
onto the observer’s null cone gives the observable quantity 
R. The luminosity distance is theoretically the same as the 
diameter distance , and is measurable provided we know 
the true absolute luminosity of the source at the time of 
emission L. If the observed apparent luminosity is l[z) then 



Let the observed number density of sources in redshift 
space be n(z) per steradian per unit redshift interval, so that 
the number observed in a given redshift interval and solid 
angle is 

ndTLdz (32) 

and over the whole sky this is 


{, 2 - - dz 

R p = mn — . 

dr 


(36) 


Thus we may substitute for R and p from ( |3l[ ) and (|36D. 

We transform ( [j^ to be in terms of redshift 2 instead 
of coordinate r by writing it as 

dftdz diR , 

-j— — -1-yy (1 -|- 2 ) -|- 47rpi?(l z) — 0 

dr dr dr^ 

and applying 

dR dz 
dz dr 
to get 


dR dz dR d R dR d z d R /dz\ 
dz dr dr ’ dz dr'^ dz'^ \drj 


dR Hz 
dz dr^ 


-b 


Hr _ , Vz 


T) (37) 


Integrating with respect to r and using (pbl) gives 


d 

dr 


dz dR 
dr dz 


V + z) 


f - dr 

dr = — 47vp{z)R{z){l-\-z) — dz{SS) 

Jo ^ 


dz dR . , 


i{z)n{z) 

R{z) 


(1 -b 2 ) d 2 


(39) 


and we used the origin conditions [{dz/dr){dR/dz)]o = 
[dR/dr)\Q = 1, and 2 ( 0 ) = 0. It follows that 

Note that this equation differs from the analogous one in 
Stoeger et al ( p3|)P| — their equation (32) — by a factor of 
(1 -b 2 ), and perhaps aptly illustrates the difference in the 
coordinate systems. To get the full model we have to solve 
the null Raychaudhuri equation (^) to get r( 2 ) (and thus 
z{r)). Equation (^) is a first integral of (§) . This has to 
be integrated one more time to obtain r(z). We must also 
specify boundary conditions at the origin r = 0, which we 
have already used in getting to (^o|): 

|(0)^J(0)f(0)^l/f(0) 

and also 


dz 

dr 


dR 

17 “ + *' 


2 ( 0 ) = 0 <b> r (2 = 0)=0. 


Airndz. (33) 

Thus the total rest mass between 2 and 2 -b dz is 


Anrhndz 


(34) 


where 7h{z) = m(T(z)) is the mass per source — i.e. the 
true density over the source number density. This primarily 
represents the evolution in the number density of sources. 
Given a local proper density p = p{t, r), and its value on the 
null cone p, the total rest mass between r and r -b dr is 


.. ‘inR^R' , 

pd‘^V = p — dr 

VTT2E 


(35) 


where d^V is the proper volume on a constant time slice, 
evaluated on the null cone. Hence by (Q, (^ and (^) 


so that, integrating dr/dz gives 
'dR, 


r{z) = 


dz 


(1 + 2 ) 


1 — 471 


i{z)n{z) 

m 


(1 -b 2 ) dz > dz . (41) 


4 THE THEOREMS 


M There they use Mq which equals fi-nrhn/R^ in the current 
notation. 


© 0000 RAS, MNRAS 000, 000-000 
















6 Nazeem Mustapha et al 


4.1 Theorem (A): 

Subject to the conditions of appendix ^ for any given 
isotropic observations i{z) & 'n{z) with any given source 
evolution functions L{z) & m{z), a set of LTB functions can 
be found to make the LTB observational relations fit the 
observations. 


4.2 Proof: — Algorithm (A): 

To obtain the LTB mass, energy and bangtime functions 
(M, E and ts respectively) from observational data and 
source evolution we would proceed as follows. 

• Take the discrete observed data for £{z, 8, (p) and 
n{z,9,(j>), average it over all angles to obtain l[z) and n{z), 
and fit it to some smooth analytic functions, such as poly¬ 
nomials. We may wish to first correct the data for known 
distortions and selection effects due to proper motions, ab¬ 
sorption, shot noise, image distortion, etc; 

• Choose evolution functions L{z) and m( 2 ) based on 
whatever theoretical arguments may be mustered; 

• Determine R{z) from L{z) and £{z) using (pl[); 

• Solve for r{z) and hence z{r), then convert func¬ 
tions of 2 : to functions of r — see appendix ^ for existence 
conditions; 

• Solve (^) and for M(r) — existence conditions 
are given in appendix 

• Determine E{r) from 

• Solve for rj from di) and dlil); 

• Solve for r(r) from and ( [l^ — L{t) and m(r) 
could now be found; 

• Determine tB^r) from (|l^). 

In practice, these equations would be solved numerically, and 
in parallel rather than sequentially, nevertheless the above 
would determine the numerical procedure within each inte¬ 
gration step. □ 

By determining the 3 arbitrary functions, we have spec¬ 
ified the LTB model that fits the given observations and 
evolution functions. This result simply asserts we can con¬ 
struct a (generally inhomogeneous) spherically symmetric 
exact solution of the field equations that will fit any given 
source observations combined with any chosen source evo¬ 
lution functions. 

We assert, without proof, that if the given observations 
and source evolution functions are reasonable, then the LTB 
arbitrary functions will generate a reasonable LTB model. 
Our definition of ‘reasonable’ is intentionally rather vague. 
By reasonable observations we obviously include the ac¬ 
tual data, suitably processed to account for selection effects. 
We also include ‘realistic’ hypothetical alternatives, but not 
functions that are wildly different from reality. Reasonable 
evolution functions are hard to define since the actual ones 
are not well known, especially at larger 2 values. By a rea¬ 
sonable LTB model, we mostly mean that the density and 
expansion rate will be within realistic ranges. A less crucial 
criterion is that there will be no shell crossings too close 
to the past null cone. Evolving the model a long time away 
from the null cone, either forwards or backwards, may intro¬ 
duce shell crossings because the data is imprecise. In general 
we don’t expect shell crossings on the large scale — i.e. two 
or more different large scale flows of galaxies in the same 


region — nevertheless it is conceivable and in that case the 
LTB description is inapplicable]^ 


4.3 Corollary (B): 

A LTB model can be found to fit the observations with 
zero evolution — m = constant, L = constant. 


4.4 Proof: 

This is an obvious consequence of (A). □ 

Given realistic data, these models will be inhomoge¬ 
neous. Indeed this is the reason that non-zero evolution 
functions have been introduced (otherwise, observations are 
incompatible with a FLRW universe). 


4.5 Theorem (C): 

Subject to the conditions of appendix for any given 
isotropic observations £{z) & n{z), and any given LTB 
model, source evolution functions L{z) & m{z) can be found 
that make the LTB observational relations fit these obser¬ 
vations. 


4.6 Proof: — Algorithm (C): 

We adapt the above algorithmic procedure to prove this. 


• As before, average the data over all angles, and fit it to 
smooth functions £{z) and n{z); 

• Specify two of the three functions M{r), E(r) and t b (r), 
the third being determined by the coordinate condition (CT. 
It seems expedient to choose M{r) and E{r). 

• Determine R{r) from the first order differential equa¬ 


tion in R and its r derivative — equation ([l2|). pi ~l| The func¬ 
tions s houl d be chosen to satisfy the origin conditions of 


section 2.2 


— see existence conditions in appendix 


• Calculate p(r) from (|^; 

• Solve for ts(r) as well as r(r) from (^, (|l^ ) and (|r^) 
with (0) defining r)(r) 


Integrate (|3ClD to get z(r) — appendix ^ gives the ex¬ 
istence conditions; 

• Use the given £{z) and n[z), to find L[z) from ( 0 ) and 
m( 2 ) from (|36|). From these and t{z) solve for L(t) and 
m{T), if needed. □ 


Again we assert that if the given observations and LTB 
model are ‘reasonable’, then the derived evolution functions 
will be ‘reasonable’. The idea is that we can vary the LTB 
model to which we fit the observations to some extent, but 
still keep the required source evolution functions within a 
‘realistic’ range. 


§§ Data can be extended through a shell crossing (^^, but not 
within the LTB formalism. 

Though we don’t strictly know the sign oi R = 
yj2M/R + 2E, it is fairly safe to assume it is positive on our 
past null cone on the large scales we are considering. 
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4.7 Corollary (D): 

Source evolution functions can be found that make the 
dust FLRW observational relations fit any observations. 

4.8 Proof: 

An obvious consequence of (C). □ 

Loosely put these theorems say 

(i) You can always fit isotropic observations with an LTB 
model, whatever the source evolution; 

(ii) If you fiddle the source evolution hard enough, you can 
fit the observations to any LTB or dust FLRW model. 

Although theorem (C) is an extreme case, and is likely 
to generate highly unphysical evolution functions if the LTB 
model is chosen arbitrarily, it is just a generalisation of (D) 
which is regularly used in an attempt to determine evolu¬ 
tion functions from cosmological observations. Theorem (C) 
highlights the dangers of this approach. 

A complication arises if the redshift is not monotoni- 
cally increasing with distance. We have seen from the well 
behaved numerical example in ( p5| ) for a parabolic case, that 
R{z) and p{z) may not be single valued, and that the R — z 
and p-z plots can loop. However, in compiling the real ob¬ 
servational data, we merely add all the galaxies we see at 
a particular redshift, to get a number count. Similarly, we 
merely take an average over the luminosities observed at a 
particular redshift, ascribing the variation to natural scatter 
in intrinsic properties and observational error, rather than 
to a multiply valued function. Thus we make R{z) and p{z) 
single valued by construction. So the data functions we are 
trying to fit may not lead to such a good model. In other 
words, assuming we succeed in constructing a well behaved 
LTB model from the data, it may not be the LTB model that 
best represents the real universe. It seems unlikely — though 
not entirely impossible — that there will be a reliable way 
of de-convolving the superposed parts of these observational 
data curves, or even of discerning whether loops are present. 
It is hard to predict how likely this scenario is. 


5 CONCLUSIONS 

We have shown that a LTB model (a Lemaitre-Tolman- 
Bondi spherically symmetric dust cosmology) can be found 
to fit any given set of observations of source counts n{z) and 
luminosity/area distance R{z), averaged over all angles, and 
any evolution functions for source luminosity L{z) and mass 
per source m{z). In other words, even if we accept isotropy, 
then demonstrating homogeneity — rather than assuming 
it must hold because of the Copernican principle — requires 
more than these observations. Conversely, our result can be 
used to determine the degree of inhomogeneity from the ob¬ 
servations and given source evolution functions. 

If the demonstration of homogeneity depends on know¬ 
ing the source evolution, and validation of source evolution 
theories depends on knowing the cosmological model is ho¬ 
mogeneous, then neither is proved. Thus we need methods of 
validating source evolution models that don’t depend on as¬ 
sumptions of homogeneity to establish the age at any given 
2 . Similarly deep cosmological distance measures that don’t 


depend on luminosity and are not influenced by source evo¬ 
lution would help pin down the cosmological model better. 
There are various promising developments, in particular: 

(a) distance measurement by Supernovae; 

(b) determinations of cosmological parameters via grav¬ 
itational lensing measurements; 

(c) accurate measurements of the Sunyaev-Zel’dovich 

effect 

(d) observation of CMB R Doppler peaks by the MAP 
and COBRAS/S AMBA [‘‘ * *| satellites. This will only deter¬ 
mine parameters in the neighbourhood of z = 1000, but is 
independent of source evolution all the same. 

(e) the increasing number of source evolution studies 
that look for tell-tale signs of early stages of galaxy evolu¬ 
tion, such as intense star formation, etc. 

Once again, the FLRW assumption is usually if not al¬ 
ways made in analyses of these effects. A re-analysis that 
permits inhomogeneity would be very worthwhile, as these 
techniques may well provide information complementary to 
the principal cosmological measures, that would help sepa¬ 
rate out the effects of cosmic evolution, spatial inhomogene¬ 
ity, and source evolution. Some of these issues are discussed 
in (||). 

In fact, it is already difficult to constrain the values of 
Ho, qo and A within a homogeneous dust model because of 
the uncertainty in source evolution, as pointed out in (§)■ 
In this case the value of A affects the time evolution of the 
scale factor, and so the deviation of the angular diameter- 
redshift relation from expectation for a A = 0 FLRW model 
could be due to non-zero A or to source evolution. Similarly 
the possible presence of non-baryonic dark matter — or for 
that matter, the possibility that gravity obeys field equa¬ 
tions other than Einstein’s — could significantly affect the 
cosmic time evolution, and introduce further uncertainty. 

The introduction of multi-colour observations does not 
resolve the problem in any simple way. If we have observa¬ 
tions in various colour bands — say U B & V — then we 
must replace the source luminosity evolution function by a 
set of evolution functions for the luminosity in each colour. 
Thus, if we find deviations of the observations from FLRW 
expectations, we still have a freedom to attribute this either 
to inhomogeneity or to source evolution. It’s true that young 
galaxies with lots of star formation are very blue. But, hav¬ 
ing introduced colour observations, and permitted evolution 
in colour, we must also admit the possibility of spatial in¬ 
homogeneities in the intrinsic colours of sources. We come 
back to the same problem — are the differences between 
observations in different colours due to source evolution or 
spatial inhomogeneity? The only difference here is that cos¬ 
mic evofution is fairly easily factored out, as the redshift is 
measured. 

We are not here asserting that the observable universe is 
inhomogeneous, nor are we suggesting that source evolution 
studies that assume homogenity are not worthwhile. 

The purpose of this paper is to emphasise that we don’t 


III Indeed a very interesting option is to use the bounds on the in¬ 
homogeneity obtained from the Sunyaev-Zel’dovich effect to con¬ 
strain the LTB model chosen, bearing in mind that this method 
still suffers from excessive error due to absorption effects. 

* * * i.e. Planck Survey 
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have unquestionable evidence for spatial homogeneity; and 
that we can’t have a good demonstration of homogeneity — 
or even homogeneity on average — without a reliable the¬ 
ory of source evolution, supported by measurements that 
are independent of cosmological model, and/or cosmic dis¬ 
tance measures that don’t depend on knowing the luminos¬ 
ity evolution of sources. Our best basis for assuming spatial 
homogeneity is the Stoeger-Maartens-Ellis theorem or “al¬ 
most EGS theorem” which says that, if the universe is 
expanding and the CMBR (cosmic microwave background 
radiation) is almost isotropic for all observers since decou¬ 
pling, then the universe is almost homogeneous, and more 
specihcally, the scale of CMBR anisotropy puts a limit on 
the degree of cosmic inhomogeneity. But this result depends 
on a weak form of the Copernican principle; and however 
convincing that principle is in general terms, we shouldn’t 
overstate it. This line of thought says that the earth is just 
another planet around the sun, but it doesn’t say all plan¬ 
ets are the same size or composition. It says that our galaxy 
and our supercluster are one among many, but allows several 
types of galaxy and considerable variety in galaxy cluster¬ 
ing. Thus the principle does not insist on uniformity on any 
scale, or even that the observable portion of the universe 
has a density particularly close to the “global average” — 
assuming we can define such a thing. And above all, while 
it may be true in the real universe, it is also possible that 
this is not so. 

We are entitled to deduce homogeneity on the basis of 
untested philosophical principles, such as a Copernican prin¬ 
ciple; but we must be quite clear what we are doing when 
we make such a deduction, and how it relates to possible ob¬ 
servational tests. This paper helps throw light on the latter 
issue. 
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APPENDIX A: A CHARACTERISATION OF 
HOMOGENEITY FROM OBSERVATIONS IN 
AN ISOTROPIC UNIVERSE 


We demonstrate the procedure in section ^ in the case where 
the input observations, after correcting for evolution, are in 
the RW form. It turns out that the LTB arbitrary func¬ 
tions assume their RW form. This amounts to a proof that 
a radially inhomogeneous dust universe is RW iff the area 
distance and number count relations as a function of redshift 
take the RW form. It is a special case of Theorem (B) where 
we assume that there is no evolution. These RW relations 
are 


^ _ 90 ^ + (1 “ '?o ) ( 1 - V2 go 2 + 1) 
Hoqo^{l + zY 

and 

‘i'lrrhn [go a -I- (1 — go )(1 — \/2 go a -I- 1)] 
3 Hogo(l-I-2)®V^2goF+T 


(Al) 


(A2) 


respectively Then we can integrate the null Raychaudhuri 
equation (|37|) once obtaining 


— — Ho{l -|- z)^y^2qoz + 1 


(A3) 


This may be integrated once again (illustrating with the case 
go < i) to obtain 


r 


1 


Ho{l 

— 250 ) 

-I- 

® In 

yi 

— 2go 

-1- 

® In 

yi 

— 2go 


^ yJ2qQZ -t- 1 
1 + 2 


/ V2go2 + 1 + \/l “ 2qo \ 

\ \/2qoWI — \/l — 2qo J 
Vi +Vl-2goy_ 


(A4) 


We continue by solving the ftrst order linear differential 
equation for M{z) (the effective gravitational mass) 0- 
This equation may be written as 

r Mjz) 
dz (1 + z)dR / dr 


2-Kmn 2-Kmn 

+ ■ 


{1 + z) (1 + z){dR / drY 


(A5) 


We substitute the RW area distance and number count func¬ 
tions into this and find that 

M{z) = HlqofV(l + zf (A6) 

and from ( p^ it follows that 

2E{z) = {l-2qo)Hlie{l + zf (A7) 


These two relations show that M oc (2E)®^^. We next show 
that this universe has a simultaneous bangtime. We need 
(Al) in addition to ( p+j ) and Restricting ourselves to 
the case go < it follows that 


710(1 — 250) 1 + 2 


y/2qoz + 1 


go 


^ Y2qQZ + 1 + \/l — 2qQ 

\/l — 250 y ^2502 + 1 — \/l — 250 


:ln 


(AS) 


and thus 


tfl(r) 


Ro(l - 2go) 

go ( 1 - yi - 2go ^ 

yi - 2go Vl + yi-2goy 


(A9) 


which means that the bang surface is simultaneous !* * ^ 
This, together with M oc (2E)®/^ is all we need to show. 


APPENDIX B: CONDITIONS FOR EXISTENCE 
OF SOLUTIONS 


B1 Existence of solutions r{z) and z(r) to equation 

( 0 ) 

We know R{z) —> 0 as 2 ; ^ 0, but from ( js^ ) we expect 
n{z) ~ fY{z), assuming p, rh{z) and dz/dr all ~ constant 
as z ^ 0, and of course m, n, R and (1 + 2 ) must all be > 0. 
However the existence condition is less stringent. 

Assume that near z = 0, mn{l + z)/R — S{z)z’’, where 
(J is a constant and 5(0) is finite and non-zero. Then, to 
leading order near 2 = 0, 


h{z)= / ■^(l+^)dz= S(z)z‘^dz 

Jo R Jo 


h = 


5(0) 


(o- + 1) 


which exists provided cr > — 1. Since all the terms in the 
integrand are positive, Jr is monotonically increasing. 

We expect R{z) to increase from 0 to a maximum, at Zm 
say, and then decrease asymptotically towards 0. We assume 
no looping, i.e. R{z) is single valued and df/dz doesn’t di¬ 
verge, and that there is only one maximum. Thus dR/dz 
goes from positive to negative values, and approaches 0 
asymptotically. It is evident from ( js^ ) that 472 / 1 ( 2 ) = 1 
where dR/dz = 0 in LTB models. Therefore we write 
dR/dz{l + 2 ) = P{z){z — ZmV 8 'iid {1 — 472 / 1 } = Q( 2)(2 — 
2m)^, where a f3 are constants, and Pm = P{zm) & 
Qm = Q{zm) are finite and non-zero. Then to leading order 
near z = Zm, (^) is 


r — r 


m 



2m)“ ^(tz 


Pm (2 - 2^)"-^+^ 
Qm (Q! — // + 1) 


and so r( 2 ) exists for a — // + 1 >0 

Our conditions for the existence of r{z) are: 

(i) rh, n, R and (1 + 2 ) are > 0, 

(ii) near 2 = 0, rhn/R ~ z” with cr > —1, 

(iii) dR/dz is finite everywhere, 

(iv) near 2 = Zm, dR/dz{l-\-z) ~ (z — ZmV 8'iid {1 — 472 / 1 } ~ 
(2 — 2m)^, with a — // + 1 > 0. 

The condition for the existence of 2 (r) is 

(v) z{r) is monotonic. 


* * * We could set t(0) = to in which case tg{r) = 0, but this is 
not necessary. 
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Conditions (i) & (ii) are manifestly reasonable. Condi¬ 
tions (ii^ & (iv) are more problematic. As was shown previ¬ 
ously (p5|), large enough inhomogeneities can create maxima 
and minima in z(r) and so make r{z) multi-valued, especially 
near dR/dz = 0, in which case neither (iii) nor (v) would be 
satified. However, a multi-valued r{z) manifested itself in a 
R{z) graph that looped. In practice, we don’t expect to get 
a looping R{z) from the observational data. The values of i 
and n at each a are averages over all measured values, and 
so are single valued by construction. Also z{r) was always 
single valued in the numerical examples considered in (|25|), 
so, if r{z) exists, then inverting it should not be a prob¬ 
lem. Unfortunately (iv) is unlikely to be satisfied exactly for 
real data — the maximum in R will not be at exactly the 
same value as the locus of dvr/i = 1, so one would have to 
tweak the fitted function to obtain a numerical solution. In 
other words, the function r{z) is sensitive to observational 
error here. This however is not too serious, since there will 
be a measure of freedom in the smooth functions i(z) and 
n(z) that are fitted to the discrete data. In fact this problem 
exists even if the universe were genuinely homogeneous — 
even if we knew the source evolution functions exactly, the 
best-fit curves obtained from imprecise observational data 
would need adjustment to obtain a solution. 


the apparent horizon R = 2M. Thus M{r) exists in the 
neighbourhood of r™. 

For completeness we consider zz yf 1, working to leading 
order. 


_ „-Fm(r-rmP ''/ (1-v) 


M = e 


G{r){r — Tm) ''^dr 


M = e 






For Q < V < 1, we again get 

m 


to leading order, which is the expected value. But for v > 1 
we get a divergence at r = r^. 

Thus our conditions for existence of M{r) are that 

(i) m, n, R and dz/dr are > 0, which ensure p > 0, and 

(ii) R(r) = R{z{r)) has a power-law maximum of the form 
7? ~ (r — rm)°‘ with 1 < a < 2, 

with a quadratic maximum being the most reasonable. 


B2 Existence of solutions M{r) to equation (14) 

Eq (0 has the form of an inhomogeneous linear hrst order 
ODE, 

— -h a(r)M = b(r) 

dr 


which has the formal solution 


M ^ fi- 


Mmyi-m + / bfidr 




where Mm = M{rm) and pm = p(rm). 

However we know that dR/dr goes through 0 at the 
maximum of R{r) — at Tm say, so both a{r) and b{r) are 
divergent there. It is evident that a &ib are hnite everywhere 
else, so we just have to show M (r) exists in the neighbour¬ 
hood of this divergence. Suppose that, near r = rm, dR/dr 
is of the form dR/dr ~ (r — rmY, where v > Q and is con¬ 
stant, so a(r) = F{r){r — rm)~^ and b{r) = G(r)(r — rm)~‘', 
where F{r) and G(r) are finite, positive and non-zero. We 
expect zz = 1. Then to leading order near Cm, for zz = 1, 




where Em = E{rm) 

, SO 

M = {r - rm)~^"^ 

pr 

0+ G(r)(r 

Fm 

M = {r - rm)~^"^ 

Gm{r - rm)^”‘ 

Fm 


and thus 

M= ^ 
Fm 


to leading order. Comparison with (^^ shows Mm = Rm/‘2, 
which is consistent with the fact that dR/dr = 0 lies on 


B3 Existence of solutions R{r) to equation (12) 

The equation is 


dR 


2M 


—— — \/1 + 2E — A / —^— -f 2E 
dr y R 

assuming that we take the positive root on the right — i.e. 
that large scale recollape has not begun anywhere on our 
past light cone. Near r = 0, E, M R aXl go to 0, but our 
origin conditions require M ~ r^, E r^ & 7J ~ r, so the 
solution exists here. Where R = 2M the r.h.s. is zero, so 
R has a maximum. We already have 2E > — 1 for a well 
behaved metric, and 2M/R & 2E are separately positive for 
parabolic and hyperbolic models, while for elliptic models 
we see from (Ei fc (0) that {—2E)R/M = (1 — cosp) < 2, 
so \/2M/R + 2E = \/ M/R \/2 + 2ER/M is always real. 

Our only conditions for R{r) to e xist are: 

(i) the origin conditions of section 2.2 are satisfied. 


B4 Existence of solutions z{r) to equation ( ^ ) 

The origin conditions ensure that, near r = 0, d?R/dr^ ~ 0, 
dR/dr ~ 1, and p ~ constant, so that the integral exists in 
this neighbourhood. 

Where the null cone crosses the apparent horizon, R = 
2M, we have dR/dr — 0. However, we find from (0) & (g) 
that the integrand of (|3C|) is 


d'^R 

dr'^ 


+ AirpR 


f m 
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V1 + 2E - xl‘^ + 2E 
R 


/ mViT 


2E 


M' 




-E' - V1 + 2S - 

R 


R 


/ ( Vl + 2E\I^ + 2E 

R 


2M 


/ VT+2E - J^ + 2E 


2M 


R 


MV1 + 2E 

_R2 R 


/ 


(^l + 2E 



which is well behaved at i? = 2M. 

Our conditions for existence of ln(l + z) are merely 
(i) the origin conditions, E ^ V, M ~ r® near r = 0. 
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